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FREE- FLIGHT MEASUREMENTS OF DYNAMIC STABILITY 
DERIVATIVES OF A BLUNTED 120® CONE IN HELIUM 
AT MACH NUMBER 15.4 

By Frederick W. Gibson and James E. Carter 
Langley Research Center 

SUMMARY 


An investigation has been made to determine the dynamic stability derivatives of a 
blunted 120® cone for various amplitudes of oscillation and center-of-gravity locations 
at a Mach number of 15.4 in helium. The data were obtained by utilizing a free- flight 
method and these data are compared with those predicted by unmodified Newtonian theory 
and experimental results obtained in air. 


INTRODUCTION 

Considerable effort has been expended (e.g., see refs. 1 to 3) to understand under 
what conditions helium data can be expected to be indicative of results that would be 
obtained in air. Nevertheless, it is not known whether hypersonic helium data can be used 
in the prediction of supersonic characteristics of large-angle blunted cones in air. The 
aim of this investigation was to examine the air-helium simulation problem for a 120® 
blunted-nose cone as well as to provide additional dynamic stability data for this configu- 
ration. A wind-tunnel free-flight technique was employed in a Mach 15.4 helium flow to 
obtain stability coefficients free of sting effects. This technique is reported in detail in 
references 4 to 7. The models were constructed of a fiber-glass shell filled with plastic 
foam and had a lead slug extending from the nose to various locations along the axis of 
symmetry. The models were launched at various angles of attack and propelled upstream 
by a pneumatic launcher. High-speed motion pictures were taken of the flight upstream 
and return, and from these the basic time -dependent attitude data were obtained. 

SYMBOLS 

Cj 5 drag coefficient, 

^ q S 


Cl 


Lift force 
q S 

Aqo 


lift coefficient, 



Cl 


a 


lift- coefficient slope per radian 


'rxi 


pltchlng-moment coefficient, PltcM-g moment 

<I«Sd 


CniQ, pitching- moment coefficient slope per radian 


Cmq + Cjn^ dynamic- stability coefficient, assumed constant over an oscillation cycle, 
fl 0d „ ad 



m 

n 


P 

Poo 

Pt,CO 

^CO 

r 


pressure coefficient, — 

^OO 

model base diameter 

model moment of inertia about a transverse axis at center of gravity 

free- stream Mach number 

mass of model 

number of cycles 

static pressure 

free- stream static pressure 
free- stream total pressure 

free- stream dynamic pressure 
radius 

radius of base 
radius of nose 
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5 model base area, nd^/4. 

t time 

V model velocity relative to medium 

Vj model velocity relative to inertial system 

free-stream velocity 

X,Z axis system fixed with respect to moving medium 

X model position relative to medium along tunnel center line 

Xgg distance from model nose to center of gravity divided by the total body length 

z model position in vertical direction relative to tunnel center line 

a model angle of attack as defined in figure 7 

y ratio of specific heats 

6 pitch angle as defined in figure 7 

9 root mean square of 6 as defined by equation (8) 

6^ value of 0 at t = 0 

0JJ value of 0 at the end of n cycles 

^oo free-stream density 

CO circular frequency 

An arrow over a symbol denotes a vector. 

A dot over a symbol denotes a derivative with respect to time. 
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APPARATUS AND TESTS 


The tests were conducted in the 24-inch-diameter (60.96-cm) nozzle of the Mach 15 
helium flow apparatus at the Langley Research Center. This blowdown-type tunnel, 
which operates at an average test-section Mach number of 15.4, is described in refer- 
ence 8. For the present study, the stagnation pressures ranged from approximately 
234 to 330 psia (1.61 to 2.28 MN/m2), stagnation temperatures from 70° to 80° F (294° to 
300° K), and Reynolds numbers from 1.65 x 10® to 2.33 x 10® per foot (5.41 X 10® to 
7.64 X 10® per meter). 

The launcher which was used to propel the models upstream is shown in the photo- 
graph of figure 1 and schematically in figure 2. The model was launched at a predeter- 
mined initial angle, velocity, and tunnel stagnation pressure. 

MODELS 

The models were blunted 120° cones with a base diameter of 3 inches (7.62 cm) 
and with nose and shoulder radii 0.125 and 0.04 times the base diameter, respectively. 

The model construction was as shown in figure 3, that is, foam-filled fiber- glass shells 
with a lead slug symmetrical with the axis of symmetry for ballast and to change the 
center- of- gravity location. The launch angle of attack was built into the model as indi- 
cated in figure 3, An analytical balance was used to determine the mass of the models, 
and the moment of inertia was determined by use of a torsion pendulum. The moment of 
inertia was calculated for calibration bodies and the periods of these bodies were deter- 
mined experimentally by use of the device shown in figure 4. The periods of the models 
were then obtained in a similar manner, and the moments of inertia were obtained from 
the calibration curve of I as a function of the square of the period for the calibration 
bodies. The positions of the centers of gravity were obtained by use of the device shown 
in figures 1 and 4. The physical properties of the models and test parameters used in the 
experiment are given in table I. 


DATA REDUCTION 

The raw data were obtained from motion pictures taken of the flight in both the ver- 
tical and horizontal planes, and consisted of time histories of translational position and 
oscillation amplitude of the model. A series of motion-picture frames from a typical 
test is shown in figure 5, and a plot of the data from another flight is shown in figure 6. 
For the tests reported herein the model motion was planar, although the plane of motion 
was inclined with respect to the horizontal and vertical planes, as shown in the typical 
photographs of figure 5, Only the tests for which at least two cycles of planar motion 
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were observed are reported herein. The data acquired should be considered qualitative 
only, as it was not possible to determine the degree of uncertainty in the measurements 
because of the small damping which the models demonstrated and the difficulties 
encountered in repeating tunnel operating conditions. 

The data- reduction technique used in this analysis is described in reference 7 and 
a brief description is given in the following sections. For convenience, a detailed deri- 
vation of the equations used to obtain the experimental drag and stability coefficients is 
given in the appendix. 


Drag 


The coordinate system used for the data reduction references the model position to 
the moving gas; figure 7 (from ref. 7) shows this coordinate system. The distance 
between the model and the moving medium is x and it is the independent variable for 
the equations of motion. The drag coefficient can be obtained from the translational 
equation of motion: 


m 


d^x 

dt2 


"ip v2sCd 

^ OO 


( 1 ) 


After rearranging equation (1) and changing the independent variable from time to dis- 
tance, the following equation for the drag coefficient results: 


Cd - 


2m 

'p S 

GO 


in 11+^ 

» nr. 


dx 


( 2 ) 


A linear fit through the data for In 1 + 


V. 


as a function of 


coefficient valid over several cycles of the oscillatory motion, 
against x is shown in figure 8. 


in 1 + 


X yields an effective drag 
A representative plot of 


Static Stability 

The solution to the linearized equation of planar angular motion (i.e., linear aero- 
dynamic coefficients Cj^ = and Cj) = Constant) is 


0 = cos 


p Sd 

OO 

21 


1/2 

Cm„ + ^ ) X 


(3) 
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where 0 is a function of the distance x and 


p S 

^OO 

4m 






In general, 


p Sd „ 

— — Cm » and hence equation (3) becomes 
21 « 


0 = cos 


p Sd 

oo 

21 


m 


a 


1/2 


( 4 ) 


If, in equation (4), x is replaced by V^t, then Cjjiq, can readily be related to the cir- 
cular frequency of motion w as follows: 


C 


niQ, 



( 5 ) 


which yields an effective constant pitching- moment slope valid over several cycles of the 
oscillatory motion. 


Dynamic Stability 

The experimental dynamic- stability coefficient Cniq + Cm^ was obtained from the 
amplitude envelope 


0 = 0Qe 


Ax 


( 6 ) 


Inserting the definition of X in equation (6) and evaluating at the end of n cycles (i.e., 


replacing x with x - 


27rnV. 


0) 


and setting 0 ~ 6^ leads to the following equation: 




+ C 


m 


a 




Trn 


1 '"n ^ 

In — + 

do 




( 7 ) 


It should be noted that values of Cl^ were not obtained in the present experimental 
investigation. However, an experimentally determined value of Cl^ was taken from 
reference 9. This value was obtained in air and is discussed subsequently. 


Data-Correlation Parameter 

The root mean square or the amplitude envelope was used to correlate the data. 
If the root mean square is defined as 0, 
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( 8 ) 




X 


Substitution of equation (6) into equation (8) and integrating through n cycles yields 


e = 


e 2 e 2 

^n - *^o 


i^o) 


1/2 


(9) 


THEORETICAL CONSIDERATIONS 
Air- Helium Simulation 

The air-helium simulation problem has been considered (e.g., see refs. 1 to 3), but 
at the present time the findings are still incomplete, particularly for large-angle blunted 
cones. Figure 9 presents the calculated pressure distribution at zero angle of attack on 
the surface of the models for helium flow at = 15.4 and for the air flow at 
Moo = 8.0, 4.0, and 5.0. These theoretical calculations were made with the computer pro- 
gram given in reference 10. The agreement between the air results at = 3.0 and 
the helium results at = 15.4 is excellent. This result partially justifies the use in 
equation (7) of the experimental value of obtained in air, and correspondingly 

allows the comparison of stability data for a 120° blunted cone under the same conditions. 

Stability- Derivative C alculations 

The theoretical dynamic and static stability derivatives were calculated by using 
equations developed from those in reference 11, which are based on unmodified Newtonian 
theory. It was found that the contribution of the spherical nose to these derivatives was 
very slight and, hence, could be ignored. The use of Newtonian theory is restricted to 
low values of oscillation frequency. In view of the relatively high frequencies of the 
present experiment, it was not expected that the dynamic stability derivative would be 
predicted accurately. 


RESULTS AND DISCUSSION 

The aerodynamic stability and drag data of the tests and the theoretical calculations 
are presented in tables I and II and figures 10 to 13. The drag data, as shown in figure 10, 
agrees well in magnitude both with that predicted by unmodified Newtonian theory and 
with the experimental result from reference 9 obtained at Mach 2.96 in air. In figures 11 
and 12, the dynamic- stability coefficient C^q + is plotted as a function of mean 

amplitude and center- of- gravity position, respectively, and compared with the Newtonian 
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results over the range of the experiments. The experimental values are considerably 
greater in magnitude than the Newtonian values; however, the present data is seen in fig- 
ure 12 to have fair qualitative agreement with that from reference 12 (air data) extrap- 
olated to the present center- of- gravity locations. 

The static stability derivative is plotted as a function of center- of- gravity position 
in figure 13. The agreement with Newtonian theory is fair; however, the scatter in the 
experimental data precludes any conclusion that the static stability decreases as the cen- 
ter of gravity moves aft in the model - a trend which is predicted by Newtonian theory. 
Shown as a solid line in figure 13 are the data obtained from reference 9 and 

transferred to the present range of center of gravity. The lack of better agreement with 
the present results is unexplained, although it is probably indicative of the breakdown of 
the air -helium simulation for a large -angle blunted cone at an angle of attack. 

CONCLUDING REMARKS 

The results of a study of the stability derivatives for a blunted 120° cone at 
M<,o = 15.4 in helium have indicated that the free-flight technique employed is a useful tool 
for obtaining data free of sting effects. However, scatter in the data, particularly for the 
dynamic- stability coefficient, was too great to permit conclusions to be drawn regarding 
trends with center- of- gravity location and amplitude of oscillation. The measured values 
of the dynamic- stability coefficient were considerably greater in magnitude than the values 
predicted by unmodified Newtonian theory. 

The lack of agreement of the present static- stability data obtained in helium with 
those previously obtained in air by other investigators (NASA TN D-4719) indicates that 
simulation of supersonic aerodynamic data with hypersonic helium data is not an effica- 
cious technique for large- angle blunted cones, 

Langley Research. Center, 

National Aeronautics and Space Administration, 

Hampton Va., July 29, 1970. 
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APPENDIX 


DERIVATION OF EQUATIONS 
Drag Coefficient 

The experimental drag coefficient is found from an equation which results from 
modifying the translational equation of motion: 

dt2 2 00 

(iV 1 

»m — = r P V oCn 
dt 2 ^ ^ 


«m 


= 1 p VSCj) 

dt 2 ^00 D 


(Al) 


and since 


then 


_ dx d -yd 
dt dt dx dx 


d(ln V) 1 _ 

„m — — — = ^ p oCt% 

dx 2 oo ^ 


or 


D 


2m d(lnv) 
"p S dx 


This equation can be rearranged as idllows: 


d(ln V) _ 

+ Vi)_ 

d|ln 

W-SJ 

/ _d(lnV«,) _ ^ 


d 


dx 

dx 


dx 

dx 

dx 


dx 


Hence, the final drag-coefficient equation for data reduction becomes 


d 

ln(l + Zl] 

2m 



p S dx 

oo 


(A2) 
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Appendix - Continued 


Dynamic -Stability Coefficient 
The equation of planar angular motion is 

I ^ - i p v^Sd (Cjn + f- - 4 P V%dC^ = 0 (A3) 

2 °° \ ”^q \V dt j 2 °° ^ 

Changing the independent variable from t to x and using the x translational equation 
(Al) gives 


d _ dx d -yd 
dt dt dx dx 

d^ _ Y d \ _ Y dV y 2 .d _ d ^ y 2 d^ 

^j2 dx \ dxj dx dx 2m dx ^2 

After assuming linear aerodynamics (that is, Cl = Cl„Q!, Cj) = Constant, and 

\ V O' 

Cm = CmQ,<3!)j the equation of motion becomes 


d^e 

dx^ 


P s 

^ 00 

2m 


Fd + 


md^ 

I 



dx 


p Sd 
00 

21 




= 0 


(A4) 


It is now necessary to substitute a = f(0); this function is easily derived from the fol- 
lowing equation which results from a small-angle assumption: 


a = 6 - 


1_ ^ 

V dt 


(A5) 


If 


gravity is ignored, the vertical translational equation is 
dt2 2m 


(A6) 


Since a and 6 are nearly equal, it is permissible to set 
Cl = Cl(P) “ T ^ ^ 


which results in 
dZ 


dt 


2m 


I 


Y^e dt 


(A7) 
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Appendix - Concluded 


Since the motion is lightly damped, an approximate equation for angular motion is obtained 
from equation (A3) by setting Cmq + = 0, Again, since a and 6 are nearly equal, 

^m the resulting equation is: 


,d2s P.V^Sd 
dt2 " 2 




a. 


(A8) 


Solving for 6 from equation (A8) and substituting into equation (A7) yields 

(A9) 


(AlO) 

If equation (AlO) is substituted into equation (A4), the equation of motion becomes the 
following second- order differential equation with constant coefficients: 


dz _ I Lq! de 

dt mdCniQjdt 

and hence equation (A5) becomes 
md Cma V dt 


d20 

^ 2m 


CD-CL„,-^^(Cm,+ C^4 


P Sd 

d6 00 p 0 — 0 


(All) 


The solution to this equation is 


B = ^ cos 


p Sd 
21 


1/2 


(A12) 


where 


p S 


X = 


4m 


Cj) - Cl_. + ~~ i C m.. + C 


cn I 


'm„ <-m 


a 
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TABLE I.- MODEL PROPERTIES AND TEST CONDITIONS 
d = 0.25 ft (0.0762 m); S = 0.049 ft^ (0.00455 m^); C-r =-1.29 (from ref. 9) 


Case 

m 

^cg 

I 

I 

Pt,0O 


Cr, 

slugs 

kg 

slug-ft^ 

kg-m^ 

md^ 

psi 

kN/m^ 

psf 

N/m2 

1 

4.23 X 10-3 

0.0617 

0.45 

5.91 X 10-6 

8.04 X 10-6 

0.0224 

262 

1806.5 

130 

6224.4 

1.55 

2 

5.42 

.0791 

.45 

7.42 

10.09 

.0219 

294 

2027.1 

146 

6990.4 

1.54 

3 

5.39 

.0787 

.50 

7.77 

10.57 

.0231 

294 

2027.1 

146 

6990.4 

1.57 

4 

5.39 

.0787 

.48 

7.22 

9.82 

.0214 

322 

2220.1 

160 

7660.8 

1.54 

5 

5.39 

.0787 

.48 

7.20 

9.79 

.0214 

313 

2158.1 

155 

7421.4 

1.56 

6 

3.36 

.0490 

.36 

5.43 

7.38 

.0259 

330 

2275.3 

164 

7852.3 

1.41 

7 

3.32 

.0485 

.36 

4.74 

6.45 

,0228 

313 

2158.1 

155 

7421.4 

1.38 

8 

3.32 

,0485 

.40 

5.52 

7.51 

.0266 

290 

1999.5 

144 

6894.7 

1.54 

9 

5.46 

.0797 

.48 

8.12 

11.04 

.0238 

234 

1613.4 

116 

5554.1 

1.76 

10 

3.27 

.0477 

.55 

5.00 

6.80 

.0245 

235 

1620.3 

117 

5602.0 

1.40 









.rn = O.I25d 


^o\N~ density plastic foam 


Fiber-glass shell 


■ Lead slug 


Built-in launch 
angle 


Shoulder rodius = 0.04d 


Figure 3.- Dimensions and construction of the models. 
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Sharp points 



Precision bearing 


Center-of~ gravity locator 


Torsion pendulum 

Figure 4.- Devices for measuring moment of inertia and locating 

center of gravity. 



L-70-4743 

Figure 5.- Photograph showing a typical flight upstream in the test section 
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Figure 10.- Experimental and calculated values of the 
drag coefficient. 
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Figure 11.- Dynamic- stability coefficient as a function of 
root-mean- square amplitude. The theoretical values 
encompass the range of center-of-gravity locations of 
the experimental data. 
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Figure 12,“ Experimental and theoretical values of the dynamic “Stability 
coefficient as a function of center“ of “gravity location. 
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Figure 13.“ Variation of experimental and theoretical values of 
static “Stability coefficient with center“ of “gravity location. 
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